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ABSTRACT 

The recent measurement by WMAP of a large electron scattering optical depth Te — 0.17 ± 0.04 is 
consistent with a simple model of reionization in which the intergalactic medium (IGM) is ionized at 
redshift z ~ 15, and remains highly ionized thereafter. Here, we show that existing measurements of 
the IGM temperature from the Lyman-alpha (Lya) forest at z ~ 2 — 4 rule out this "vanilla" model. 
Under reasonable assumptions about the ionizing spectrum, as long as the universe is reionized before 
z = 10 and remains highly ionized thereafter^ the IGM reaches an asymptotic thermal state which is 
too cold compared to observations. To simultaneously satisfy the CMB and Lya forest constraints, the 
reionization history must be complex: reionization begins early at z > 15, but there must have been 
significant (order unity) changes in fractions of neutral hydrogen and/or helium at 6 < z < 10, and/or 
singly ionized helium at 4 < z < 10. We describe a physically motivated reionization model that satisfies 
all current observations. We also explore the impact of a stochastic reionization history and show that 
a late epoch of (Hell— i-Helll) reionization induces a significant scatter in the IGM temperature, but the 
scatter diminishes with time quickly. Finally, we provide an analytic formula for the thermal asymptote, 
and discuss possible additional heating mechanisms that might evade our constraints. 

Subject headings: cosmology: theory — intergalactic medium — quasars: absorption lines — cosmic 
microwave background 



1. INTRODUCTION 

The detection by the Wilkinson Microwave Anisotropy 
Probe ( WMAP) of a large optical depth Te to electron scat- 
tering has opened a new window in studies of the ultra- 
high redshift (z ^ 15) universe (Kogut et al. 2003, Spergel 
at al. 2003). Taking the value = 0.17 ± 0.04, inferred 
from the polarization and temperature anisotropics of the 
cosmic microwave background (CMB), at face value, im- 
plies the reionization of the universe must begin very early. 

The optical depth to electron scattering is given by (e.g. 
Dodelson 2003) 
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with hydrogen, the reionization redshift changes slightly 
to z = 15.3 ±2.7. 1 

Two features are noteworthy. First, helium reionization 
has a sub-dominant effect on t^. compared to hydrogen. 
Second, since the electron scattering optical depth is con- 
trolled by the free electron density, it is insensitive to the 
neutral fractions of hydrogen {Xm = 1 — Xhii) and helium 
(XhoI = 1 — AThcII — Xnciu), as long as they are small. 

In contrast, the (hydrogen) Lyman-alpha (Lya) optical 
depth, inferred from the spectra of distant quasars, is ex- 
tremely sensitive to small amounts of neutral hydrogen. 
The Lya optical depth at mean density is (Gunn & Peter- 
son 1965) 



where H(z) is the Hubble parameter at redshift z, ax is 
the Thompson cross-section, and ne(z) is the proper free 
electron density. This can be rewritten as 
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where fif,, and f2„j are the baryon and matter densities in 
fraction of the critical, h is the Hubble constant in units of 
100 km/s/Mpc, and Yp = 0.244±0.002 is the helium mass 
fraction (Buries, NoUett & Turner 2001). The fractions of 
ionized hydrogen Xaii, singly ionized helium XhcII, and 
doubly ionized helium XhoIii, are functions of z. 

Ignoring helium, the observed Te = 0.17 ±0.04 is consis- 
tent with a universe in which Xhii changes from essentially 
zero to close to unity at z = 17 ± 3, and Xmi ^ 1 since 
(Kogut et al. 2003). If helium is fully ionized together 



Whether Xui is 10^^ or 10~^, for instance, makes a big 
difference to Tq . Using a model for large scale structure in 
the intergalactic medium (IGM), the observed mean Lya 
transmission at z 6 (or the lack thereof i.e. the Gunn- 
Peterson trough; Becker et al. 2001) implies a (la) lower 
limit on the hydrogen neutral fraction for a fluid element at 
mean density: Xhi > 10""*. The analog for Ly/3 provides 
a stronger constraint, due to the smaller absorption cross- 
section: Xm > 5 X 10""' (taken from Lidz et al. 2002; see 
also Cen & McDonald 2002, Fan et al. 2002). ^ 

^Throughout this paper, we adopt Ui,h'^ = 0.024, Qmh^ = 0.14, 
and h = 0.72 the central best-fit values measured by WMAP 
(Spergel et al. 2003). 

^Note that some authors quote volume— or mass-weighted neutral 
fractions. We here use the neutral fraction for a fluid element that 
happens to be at the cosmic mean density, motivated by the fact 
that we will discuss the temperature evolution of fluid elements with 
the same property in this paper. 
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Taken at face value, these two separate observations are 
therefore consistent with the simplest "vanilla" model in 
which the universe is reionized in a single step: the neutral 
fraction experiences a drop of order unity at at 2; > 15, and 
it remains <C 1 thereafter. Our goal in this paper is to con- 
front this model with other existing observations, and to 
see if additional constraints can be put on the reionization 
history of our universe. 

Several authors have pointed out that the evolution of 
the Lya optical depth suggests reionization might take 
place not much earlier than z ^ 6 (Becker et al. 2001, 
Djorgovski 2001, Gen & McDonald 2002, Gnedin 2001, Ra- 
zoumov et al. 2002; but see also Barkana 2001, Songaila 
& Cowie 2002). This is based on an extrapolation of 
the mean transmission measurements from lower redshifts. 
The small number of lines of sight used (a single quasar 
was employed for the measurement at 2: ~ 6; but see Fan et 
al. 2003 for three new z > 6 sources with Gunn-Peterson 
troughs), the challenge of sky subtraction and continuum 
extrapolation in Gmm-Peterson trough measurements (see 
discussion in Becker et al. 2001; see also Hui et al. 2001), 
as well as our still maturing understanding of radiative 
transfer during reionization, motivates us to look for other 
clues for a late period of reionization. 

The main idea is quite simple. Reionization typically 
heats up the IGM to tens of thousands of degrees, and the 
gas subsequently cools due to the expansion of the uni- 
verse as well as due to other processes. If the universe was 
reionized early, and has stayed highly ionized thereafter, 
photo-ionization heating of the gas cannot overcome the 
overall cooling, and the IGM might reach too low a tem- 
perature at low redshifts compared to observations. This 
idea is not new (e.g. Miralda-Escude & Roes 1994, Hui 
& Gnedin 1997, Haehnelt & Steinmetz 1998, Hui 2000, 
Theuns et al. 2002). Theuns et al. (2002), in particular, 
deduced a limit on the reionization redshift {z < 9) based 
on temperature measurements by Schaye et al. (2000), 
assuming a quasar-like ionizing spectrum, and that Hell 
reionization occurs at z ~ 3. Our objective here is to seek a 
formulation of this argument that is as clean and robust as 
possible, that reveals clearly the underlying assumptions, 
and to check the consistency with the IGM temperatures 
of specific models that produce the high value of Tg mea- 
sured by WMAP. 

The Lya forest temperature measurements we will use 
are taken from Zaldarriaga, Hui & Tegmark (2001; ZHTOl 
thereafter): Tq = 2.1±0.9xl0'' K at z = 2.4, 2.3±0.7xl0'' 
K at ^ = 3, and 2.2 ± 0.4 x 10'' K at z = 3.9. Note that 
2 (J errorbars are quoted, and the temperature Tq was de- 
rived for fluid elements at the mean density, consistent 
with our modeling of the IGM temperature in the rest 
of this paper. There have been several other measure- 
ments in the past (Ricotti, Gnedin & Shull 2000, Schaye 
et al. 2000 [STOO], Bryan & Machacek 2000, McDonald et 
al. [MMOl], Meiksin, Bryan & Machacek 2001). The ones 
that are easiest to compare, because they are based on very 
similar datasets, are STOO, MMOl and ZHTOl. The for- 
mer two are based on line width measurements, while the 
last one makes use of the small scale transmission power 
spectrum. A virtue of the last method is that the tem- 
perature constraints come from marginalizing over a wide 
array of parameters, including the slope and amplitude 
of the primordial power spectrum (e.g. Hui & Rutledge 



1999), and the equation of state index. It is reassuring that 
MMOl and ZHTOl, using very different methods and em- 
ploying different simulations (MMOl using hydrodynamic 
simulations, and ZHTOl using N-body simulations with a 
marginalized smoothing to mimic Jeans smoothing) , agree 
well with each other. The results of STOO are somewhat 
discrepant from these two works - the reader is referred to 
ZHTOl for further discussions. 

An important issue in determining the temperature of 
the IGM is the second ionization of helium (Hell^Helll). 
As we will see, this can be an important source of heat- 
ing at low redshifts z ~ 3 — 4. There are several lines of 
evidence that suggest Hell might be reionized at 2: 3, 
including observations of Hell patches that do not seem to 
correlate with HI absorption (Reimcrs et al. 1997, Ander- 
son et al. 1999, Heap et al. 2000, Kriss et al. 2001, Jakob- 
sen et al. 2003), increase in IGM temperature (Schaye 
et al. 2000; Theuns et al. 2002), the evolution of the 
hardness of the ionizing background spectrum (Songaila 
1998), and evolution of the moan transmission (Bernardi 
et al. 2003). On the other hand, it is unclear if the fluctu- 
ations in Hell absorption observed in a few lines of sight 
might not be due to a naturally fluctuating IGM (Miralda- 
Escude, Haehnelt & Rees 2000); the IGM temperature 
measurements by MMOO and ZHTOl are consistent with 
no feature at z ~ 3; power spectrum evolution also seems 
to argue against Hell reionization at 2 ~ 3 (McDonald & 
Scljak, private communication). In this paper, when we 
consider constraints from the temperature measurement, 
primarily at z ~ 3.9, we therefore leave two options open; 
Hell can be ionized or not ionized by 2 = 3.9. With this 
explained, we can now state our vanilla model in more 
concrete terms. It has two variants: 1. both hydrogen 
and helium are fully reionized at z > 15, and they remain 
highly ionized thereafter; 2. the same as 1. except that 
helium is only singly ionized, and remains so until at least 
past z = 3.9. 

The paper is organized as follows. In §2, we explain the 
idea of a thermal asymptote for the IGM, and use it to de- 
rive a constraint on the hardness of the ionizing spectrum 
if the universe were to reionize before z = 10, and remains 
highly ionized thereafter. In §3, we discuss limits on the 
hardness of the ionizing spectra, and show that these spec- 
tra fall short of making the thermal asymptote sufficiently 
hot to match observations. We then work out illustrative 
examples in §4 of how order unity changes in the neutral 
fractions at z < 10 can reproduce the temperature mea- 
surements, while being consistent with the WMAP data. 
We go on to offer a physically motivated model in §5, and 
discuss the implications of the fact that different fluid ele- 
ments are reionized at different times (stochastic reioniza- 
tion). We conclude in §6. 

2. THERMAL ASYMPTOTICS 

Several thermal processes are at work in a photoionized 
IGM. They are described in detafl in Hui & Gnedin (1997). 
Here is a brief qualitative summary: 

• Adiabatic heating/cooling. Gas elements can heat 
up or cool simply due to adiabatic contraction or 
expansion. The overall expansion of the universe 
drives a temperature fall-off as (1+z)^ as z decreases. 

• Photoionization heating. Photons inject energy into 
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the gas by ionizing hydrogen or hehum. 

• Recombination coohng. Protons and electrons (or 
ionized hehum and electrons) can cool by recombin- 
ing and radiating energy away. 

• Compton cooling. At sufficiently high redshifts, z > 
10, Compton scattering of free electrons with the 
lower temperature CMB can be an important source 
of cooling. 

Typically, photoionization heating provides the domi- 
nant source of heating for the tenuous IGM. This oc- 
curs in primarily two forms. One is during what we 
call reionization 'events' - these are (possibly extended) 
periods in which a given gas element experiences order 
unity changes in the fractions of HI/Hel/Hell. The other 
is photo-heating for an already highly ionized plasma - 
this occurs through photoionization of small amounts of 
HI/Hel/Hell, whose abundances are determined by pho- 
toionization equilibrium. The former provides a big boost 
to the temperature, while the latter determines the asymp- 
totic thermal state of the IGM. Typical temperature evo- 
lutions are illustrated in several toy models in Fig. 1. 

Fig 1 shows the evolution in temperature Tq for a fluid 
element at mean density, with a variety of reionization his- 
tories. These histories have different reionization redshifts 
and initial temperature jumps, selected here only for illus- 
tration. The code for computing the thermal and chemical 
evolution is described in Hui & Gnedin (1997). For each 
thermal curve, one can see a rise to high temperatures after 
a reionization event. Thereafter, the varieties of reioniza- 
tion histories and initial reheat temperatures all result in a 
rather similar late time thermal asymptote (black dotted 
line). In fact, if reionization occurs before z ~ 10 (and 
the IGM stays highly ionized since, or to use our earlier 
terminology: no reionization events take place thereafter), 
the temperature is always very close to the asymptote for 
z < 4, which is where temperature measurements exist 
(points with error-bars). This is true even if reionization 
is not a single step process, as illustrated by the dashed 
(magenta) line, as long as complexities in the reionization 
process occur before z ~ 10. The term 'complexities' here 
has a rigorous meaning: it refers to two or more episodes 
during which the neutral fractions of hydrogen and/or he- 
lium undergo changes of order unity. In all thermal curves 
shown in Fig. 1, with one exception, such changes take 
place before about z ~ 10. The exception, i.e. the left- 
most curve where reionization takes place at z ^ 7^ is also 
the only one where the asymptote is not reached by z = 4 
- this is because there is not sufficient time for the IGM 
to cool after a recent episode of reionization. 

Fig. 1 therefore illustrates two very important general 
facts: 

1. The thermal state at z < A does not remember the 
part of reionization history prior to z ^ 10. 

2. The IGM does, however, retain short-term memory 
of reionization events in its recent past. 

What determines the late time thermal asymptote? It is 
the combination of photoionization heating, recombination 
cooling, and adiabatic cooling. An approximate analytic 



5x10" 



4x10" 



,3x10" 



o 



2x10" 




Fig. 1. — The figure shows the thermal asymptote (black 
dotted line; eq. [4]), and illustrates the fact that a wide 
range of ionization histories result in the same IGM tem- 
perature by redshift z = 4, unless reionization occurs late. 
Each colored solid line describes the evolution of the tem- 
perature (To) for a fluid element at mean density, accord- 
ing to a different reionization history (and a different ini- 
tial reheat-temperature). The (magenta) dashed line il- 
lustrates the thermal evolution for a complex reionization 
history - such complexities do not stop the temperature 
from reaching the late time asymptote, as long as they 
take place early, before z ~ 10. The points with (2 a) 
error-bars on the left are measurements of Tq from the 
Lya forest (Zaldarriaga et al. 2001). 



expression can be derived for the thermal evolution of the 
IGM (Hui & Gnedin 1997). We find that the following 
simple formula provides a more accurate (~ 5%) fit to 
the numerically computed thermal asymptote (for Tq at 
z = 2 — 4) , for a variety of spectral shapes we have tested: 



To = [B{1 + z)"-^] 
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Here, fcs is the Boltzmann constant, the symbol Ji, de- 
notes the ionizing intensity as a function of frequency (it 
has unit of energy per frequency per time per area per 
ster-radian) , ai denotes the ionization cross-section for the 
respective species {i — HI, Hel or Hell), Vi is the ioniz- 
ing threshold frequency, and h is the Planck constant (in 
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distinction from the Hubble constant h). It is useful to 
remember h^-Hi/fes = 1.57807 x 10^ K, ^Hei = 1.808 
and uueii = ^vm- 

The asymptote given in eq. (4) assumes that at late 
times, hydrogen is highly ionized. That is why the term 
due to photo- heating of HI, Thi, has the form shown in 
eq. (5). The denominator of Thi, when multiplied by 
4-77, is the photoionization rate of neutral hydrogen - its 
inverse is hence proportional to the (small) hydrogen neu- 
tral fraction under ionization equilibrium. Its numerator 
gives the photo-heating rate per neutral atom. The com- 
bination of factors in Thi therefore gives a temperature 
scale whose amplitude is proportional to the net amount 
of HI photo-heating. Similarly, the terms XhcIiThci/16, 
and 5.6XhcIii7hcIi/16 in eq. (5) quantify the importance 
of photo-heating of Hel and Hell respectively. The factors 
of Xhcii and Xneiii arise due to photoionization equilib- 
rium. If, asymptotically, helium is doubly (singly) ionized, 
then Xhoii (-^^hoIii) can be set to zero. 

It is important to emphasize that the thermal asymptote 
given in eq. (4) is determined completely by the shape (or 
'hardness') of the ionizing spectrum Jy, but not its ampli- 
tude. For a power-law J^, oc the thermal asymptote 
can be recast simply as: 

To = 2.49x lO^K X (2 + a)-T^ (6) 

assuming, asymptotically, helium is doubly ionized. 

Given the above, we can ask the following question: how 
hard does the ionizing spectrum have to be for the thermal 
asymptote to match the observed temperatures at low red- 
shifts? It suffices to use the measurement at the highest 
redshift, z = 3.9: Tq = 2.2 ± 0.4 x lO'' K (2 ct error-bar; 
ZHTOl). More concretely, what constraint can be placed 
on a if the thermal asymptote were to reach Tq > 1.6 x 10^ 
K (3 (T lower limit), by z = 3.9? It is straightforward to 
show that this requires a < 0.12, from eq. (6). 

One can also derive a similar limit if it is assumed helium 
is only singly ionized by ^; = 3.9 (i.e. XHei,-^Heiii 1, 
-^Hoii ~ 1). Assuming again a power- law Jj^ oc but 
with a cut-off for v > v^ieih the requirement is a < —2.2. 
In other words, the spectrum needs to be even harder with 
no Hell photo-heating, and a spectrum as hard as a < 
—2.2 is unrealistic in comparison with stellar and quasar 
spectra. 

To summarize: if reionization takes place at z > 10, and 
the fractions of HI/Hel/Hell experience no significant ( or- 
der unity) change after z = 1Q, a sufficiently hard ionizing 
spectrum is necessary to keep the IGM temperature high 
enough to m,atch observations at z = 3.9. Parameterizing 
the ionizing background by a power-law Ji, oc i'"", this re- 
quires a < 0.12 if Hell is reionized by ^; = 3.9, or a < —2.2 
if Hell is not reionized by then (the latter assumes J,^ is cut 
off for 1/ > t'Hoii). One should keep in mind that the only 
relevant slope of the spectrum is the slope just above each 
ionization threshold (unless the spectrum is very hard), 
because the ionization cross section cr, ~ In other 

words, a can deviate greatly from the values given above 
as long as the deviation takes place at frequencies far away 
from the ionization thresholds. Note also that the spec- 
trum described here refers to the asymptotic spectrum at 
z = 3.9. If the spectrum changes significantly after z ~ 10, 



one can view the above limits as applicable to the hardest 
spectrum betwcxai z = 3.9 — 10. 

The power index we find, a < 0.12, or a < —2.2, repre- 
sents a very hard spectrum. We next turn to the question: 
how hard can a realistic ionizing spectrum be? 

3. THE IONIZING SPECTRUM 

Two kinds of ionizing spectrum are generally discussed 
in the literature. One is quasar-like and the other is star- 
hke. 

Zheng et al. (1998) finds a quasar spectral shape of 
~ for the relevant ionizing frequencies in high reso- 

lution HST spectra. In this paper, we will follow Haardt 
& Madau (1996) and consider, conservatively, a quasar 
spectrum that goes as One should keep in mind 

that, at the relevant redshifts > 4, the known populations 
of quasars probably cannot contribute significantly to the 
(hydrogen) ionizing background. Here, we take the con- 
servative view that there might be a population of dim 
quasars that still contribute significantly to a hard spec- 
trum (Haiman & Loeb 1997). 

Stellar spectra are generally softer than a typical quasar 
spectrum. An exception is the spectrum of metal-free stars 
(Tumhnson & ShuU 2000, Bromm, Kudritzki & Loeb 2001, 
Schaerer 2002, Venkatesan, Tumlinson & Shull 2003). In 
their theoretical models, Bromm et al. (2001) find a spec- 
trum that has the following form: ^ v for v just above 
t^Hi, Ji^ ^ z/" at z/ > i^Hei, and Ji, ~ for v > uucU- 

Such a spectrum is harder than the quasar spectrum at fre- 
quencies below the Hell threshold. Note that metal free 
stars probably cause reionization early on, but it is un- 
likely they contribute significantly to the asymptotic ion- 
izing spectrum at low redshifts (Haiman & Holder 2003). 
We consider a metal free stellar spectrum here for the sake 
of being conservative i.e. assume a spectrum that is as 
hard as possible. 

The actual ionizing spectrum seen by a fluid element is 
different from the above, due to processing by the IGM. 
Haardt & Madau (1996) have done a careful calculation of 
such effects. Absorption generally hardens the spectrum^ 
However, diffuse recombination radiation from the ab- 
sorbing medium tends to compensate for this hardening. 
Haardt & Madau (1996) found that the spectrum above 
i/Heii hardens by 1.5 (to be precise, a — > a — 1.5) (see also 
Zuo & Phinney 1993). The spectrum just above uui and 
t^Hci maintains roughly the same slope as the source. 

We are therefore led to consider the following IGM mod- 
ified spectra. Let us use the symbols am, ctuei and aneii 
to denote the spectral slopes (more precisely, its negative), 
above the three relevant ionizing thresholds. A quasar-like 
processed spectrum has ani = ccnei = 1-5, and aneii = 
0.0. A metal- free-stellar spectrum gives rise to am = — 1, 
Q!Hei = 0, and aHeii = 3. At the moment (s) of reioniza- 
tion, the relevant spectrum can, however, be even harder 
(Abel & Haehnelt 1999). In principle, the spectrum can 
be hardened relative to the source by as much as a power- 
law index of 3, the 3 coming from the scaling crj a 
(see arguments in Abel & Haehnelt (1999), and also Zuo & 

•^Note the apparent paradox: absorption, by taking away ionizing 
photons, increases the photo-heating rate. This is a result of photo- 
ionization equilibrium, which makes the amplitude of the ionizing 
background irrelevant (eq. [5]). Rather it is the spectral shape that 
is important. 
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Phinney 1993). We therefore aUow the quasar-spectrum to 
have ohi = and = «HoII = —1-5, and the stellar-spectrum 
to have am — —4, and — — 3,aHcii = 1-5, at the initial 
moment (s) of reionization. This is not relevant for the 
thermal asymptote, but is relevant for the magnitude of 
temperature boosts during reionization events. 

Comparing the above spectral slopes against the limits 
obtained in the last section suggests neither quasars, nor 
metal-free stars can match the observed temperature at 
z = 3.9, if reionization occurs at z > 10, and no reioniza- 
tion event takes place afterwards. However, those limits 
were based on a strict power-law spectrum. The ther- 
mal asymptotes (eq. [4]) for our more realistic spectra for 
quasars and metal-free stars are shown in Fig. 2 and Fig. 
3 (black dotted hues), respectively. Indeed, in both cases, 
the thermal asymptotes fall short of the observed (3 cr) 
lower limit of Tq > 1.6 x 10* K at z — 3.9, regardless of 
whether or not helium is doubly ionized. 

Hence, if the universe is reionized before z ~ 10, and 
remains highly ionized thereafter, neither reasonably hard 
spectra can reproduce the IGM temperatures inferred from 
the Lya forest. 

4. ILLUSTRATIVE EXAMPLES 

The conclusion from the last section, together with 
the large electron scattering optical depth measured by 
WMAP, therefore implies that one or more of the fractions 
-'^Hi ) -'^Hci , -'i^Heii must change by order unity at z < 10, 
to give the IGM temperature boosts above the thermal 
asymptote. We therefore rule out the vanilla reionization 
models laid out in §i. One should keep in mind, however, 
this is predicated upon (reasonable) assumptions about 
the hardness of the ionizing spectrum. In this section, 
we work out some illustrative examples of what it takes to 
match the IGM temperature constraints, using the spectra 
from the last section. 

Fig. 2 shows the thermal evolution for a fluid element at 
mean density subject to a quasar-like ionizing spectrum, 
and experiencing a variety of ionization histories, for both 
the case of having helium doubly ionized (upper panel), 
and the case of only having singly ionized helium (lower 
panel; except for model E, see below). Models A and C 
both describe early reionization before z = 10, and no sig- 
nificant (order unity) changes in Xm, XhcI, ^HoII there- 
after. As explained before, both converge to their respec- 
tive thermal asymptotes by z 4, which fall short of the 
observed temperatures, especially for model C, which has 
no Hell reionization, and therefore lower temperatures. 
Model B demonstrates how early reionization can be made 
consistent with the forest temperature measurements. It 
has a late second episode of (both hydrogen and helium I 
and II) reionization at z ^ 7, which boosts the tempera- 
ture significantly above the thermal asymptote. Model E 
is similar, except that in the second episode, only Hell is 
reionized (HI and Hel were already highly ionized before 
then). It shows that Hell reionization alone can provide 
a significant boost to the IGM temperature. 

Model D is intriguing, as it shows a case where a second 
episode of hydrogen (and Hel, but not Hell) reionization 
occurs at the smallest redshift (z = 6) allowed by SDSS 
observations (Fan et al. 2002). Even with such a late 
reionization epoch, the temperature can barely match the 
observations (the temperature at z = 3.9 is 1.5 x 10~*, 
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Fig. 2. — Thermal evolution for a quasar-like ionizing 
spectrum. Upper panel: helium is doubly ionized together 
with hydrogen. The (black) dotted line shows the thermal 
asymptote (eq. [4]). The (colored) solid lines show the 
thermal evolution for two different reionization histories. 
Model A has a single-episode reionization at z = 14, and 
stays highly ionized thereafter. Model B also has early 
reionization, but experiences a drop in ionizing flux there- 
after, and starts recombining until z ~ 7 at which time 
it undergoes a second period of hydrogen and helium (I 
and II) reionization. Model B, but not A, reaches a high 
enough Tq to match observations (points with error-bars; 
ZHTOl). Lower panel: helium is singly ionized together 
with hydrogen; the double ionization of helium is never 
reached ( except in model E) . The (black) dotted line shows 
the asymptote in such a case. Models C, D (colored solid 
lines) are analogs of A and B above - the only difference 
is that here, helium remains only singly ionized. In model 
E (magenta dashed line), hydrogen and helium is (singly) 
ionized at z 17, and then Hell is ionized at z ^ 6.5. 



which is just below the 3 a limit). This is due to the 
lack of Hell reionization in this model. Therefore, with a 
quasar spectrum like the one assumed here, some amount 
of Hell reionization is necessary prior to z ^ A, to heat the 
IGM sufficiently. If the evidence that suggests that Hell 
reionization is occurring at z ~ 3 holds up (see §2), then 
this implies either that a spectrum harder than assumed 
here for quasars exists prior to z ~ 4, or else Hell reioniza- 
tion must last for a an extended period, from z 3 back 
to at least z > 4. 

Fig. 3 shows a similar exercise for a metal free stellar 
spectrum (as defined in §3). We emphasize again, how- 
ever, that it is unlikely that metal free stellar spectra re- 
main a dominant contribution to the ionizing background 
down to low redshifts (Haiman & Holder 2003). This fig- 
ure should only be viewed as an illustration of possibil- 
ities. The curves are exact analogs of those in Fig. 2. 
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Fig. 3. — Thermal evolution using a spectrum motivated 
by metal-free stars. All reionization histories are exact 
analogs of those in Fig. 2. All symbols inherit the same 
meanings. 

Two features are different from the previous figure. Model 
E, where a late period of Hell reionization occurs around 
z ~ 6, can no longer heat up the IGM sufficiently to be 
consistent with observations. This is because of the soft- 
ness of a stellar spectrum above the Hell threshold. On the 
other hand, model D, which has a late hydrogen (and Hel, 
but no Hell) reionization at z ~ 6, has no problem match- 
ing the observed temperatures, unlike its quasar analog 
shown in Fig. 2. This is because the stellar spectrum we 
adopted is, in fact, harder than the quasar spectrum for 
frequencies just above vyh and i^hoI- 

Finally, before we move on to the next section, it is inter- 
esting to explore a situation that is intermediate between 
the two extremes we have considered above i.e. Hell is par- 
tially ionized early on, and stays partially ionized there- 
after. Because the recombination rate of Helll to Hell is 
high, one might wonder if the continuous photoionization 
and recombination could give rise to a sufficiently high 
temperature."* We show in Fig. 4 the thermal evolution 
(upper panel) for a fluid element which experiences reion- 
ization at z > 10, with the Hell fraction (lower panel) 
maintaining at the 80% level thereafter. The thermal evo- 
lution assumes a quasar spectrum identical to that used 
in Fig. 2 (in particular, ckHcii = for the processed, 
post-reionization spectrum; see §3), except that a break in 
the ionizing flux just above the HeH threshold is tuned to 
give rise to Xneii ~ 0.8. Notice how the analytic asymp- 
tote (dotted line; from eq. [4], with Xyi^ii = 0.8, and 
-^Heiii = 0.2) remains an excellent approximation to the 
late time thermal evolution. One can see that the late time 
temperature is intermediate between the cases of full Hell 

*We thank Andrea Ferrara for posing the question. 



Fig. 4. — The thermal evolution (solid line, upper panel) 
using the same quasar-like spectrum as is assumed in Fig. 
2, except that the amplitude of the ionizing radiation just 
above the Hell threshold is tuned to give rise to partial 
Hell ionization (lower panel). The dotted line shows the 
thermal asymptote given in eq. (4). 

reionization and of no Hell reionization, depicted in Fig. 
2. In other words, partial Hell reionization is nicely brack- 
eted by the cases we have considered already. A caveat one 
should keep in mind, however, is that when AThoII is close 
to unity, the processed ionizing spectrum is expected to 
be still harder above the HeH threshold compared to what 
we have already assumed. Trying ancii — —1.5 (which is 
the hardest possible; see §3) in place of q;hcII = in Fig. 
4, results in a temperature at z ~ 4 of 1.3 x 10"* K, which 
still falls short of the observed temperature. For smaller 
A'hcII, one expects the spectrum to be softer. 

5. STOCHASTIC REIONIZATION HISTORY - RESULTS 
FROM A PHYSICALLY MOTIVATED MODEL 

In the above sections we have used toy models of reion- 
ization to illuminate the key issues that determine the tem- 
perature evolution of the IGM. It is interesting to con- 
sider this evolution in a physically motivated model of 
the reionization history that appears to fit all the rele- 
vant observations (including the electron scattering op- 
tical depth Te = 0.17 measured by WMAP). Based on 
assumptions about the nature and efficiency of the ioniz- 
ing sources, the reionization history can be predicted from 
"first principles" in numerical simulations (e.g .Gnedin & 
Ostriker 1997; Nakamoto, Umemura, & Susa 2001; Gnedin 
2001; Razoumov et al. 2002), and in semi-analytical mod- 
els (e.g. Shapiro, Giroux & Babul 1994; Tegmark et al. 
1994; Haiman & Loeb 1997, 1998; Valageas & Silk 1999; 
Wyithe & Loeb 2003; Gen 2003). Here we consider a semi- 
analytical model adopted from Haiman & Holder (2003), 
but modified to include the ionization of Hell^Helll. For 
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technical details, the reader is referred to that paper. In- 
clusion of Helll is conceptually straightforward. 

In this model, we follow the volume filling fractions xhii 
and XHeiii of HII and Helll, assuming that discrete ion- 
ized Stromgren spheres are being driven into the IGM by 
ionizing sources located in dark matter halos. As can be 
seen from eq. (5), photoionization of Hel (to Hell) plays 
a sub-dominant role in heating the IGM. For simplicity, 
we assume that Hel is reionized at the same time as HI, 
so there is no need to separately keep track of the filling 
fraction of Hel. Note that before the Stromgren spheres 
percolate, the radiation background is extremely inhomo- 
geneous: fluid elements inside ionized regions see the flux 
of a single (or a cluster of a few) sources, whereas fluid ele- 
ments in the still neutral regions sec zero flux."'' In this pic- 
ture, each fluid element is engulfed by an ionization front 
at a different time. In effect, each fluid element therefore 
has a different reionization history. Rather than consider- 
ing a single temperature for a fluid element at the mean 
density, it is more appropriate to consider an ensemble of 
fluid elements at the mean density, each with a different 
reionization history and temperature evolution. Note that 
this stochasticity is in addition to the IGM having a dis- 
tribution of temperatures due to density variations (see §6 
for a discussion of the latter). 

The evolutions of a^Hii and XHeiii in our model are shown 
by the solid and dashed curves in the top panel of Fig- 
ure 5 (ignore the QSO curve to the left for the moment). 
Reionization has an interesting history that reflects contri- 
butions from three distinct types of ionizing sources, and 
two different feedback effects (all of which have physical 
motivations as described in detail in Haiman & Holder 
2003). In short, ionizing sources (assumed to be mas- 
sive metal-free stars) first appear inside gas that cools 
via H2 lines, and collects in the earliest non-linear ha- 
los with virial temperature of 100 K <T< lO^K. These 
sources ionize ~ 50% of the volume in hydrogen, and 
they have sufficiently hard spectra (sec discussion above) 
that they reionize ~ 10% of the helium. However, at this 
stage (redshift 2; ~ 17) the entire population of these first 
generation sources effectively shuts off due to global H2 
photodissociation by the cosmic soft UV background they 
had built up. Soon more massive halos, with virial temper- 
atures of 10'' K < r < 2 X 10'''K form, which do not rely on 
H2 to cool their gas (they cool via neutral H excitations), 
and new ionizing sources turn on in these halos. These are 
assumed to be "normal" stars, since the gas had already 
been enriched by heavy elements from the first generation. 
These sources contimie ionizing hydrogen, but since they 
produce little flux above the Hell edge, helimn starts re- 
combining. This second generation population is also self 
limiting: gas infall to these relative shallow potential wells 
is prohibited inside regions that had already been ionized 
and photo heated to lO'^K. As a result, hydrogen reioniza- 
tion starts slowing down around z ~ 10 (see the solid curve 
in the bottom panel of Fig. 5). However, at this stage, still 
larger halos with virial temperatures of T > 2 x lO^K start 
forming. These relatively massive, third generation halos 

^These statements would no longer hold if the early ionizing 
sources had a hard spectrum extending to >1 keV energies; an inter- 
esting possibility (e.g. Oh 2001, Venkatesan, Giroux & Shull 2001) 
that we do not consider in this paper. 




Fig. 5. — Top panel. The evolution of the volume filling 
fractions Xhu and a;Heiii of ionized hydrogen and helium 
in a physically motivated reionization model. The solid 
curve corresponds to HII, and the dashed curve to Helll 
regions. The extra dashed curve to the left (QSO) corre- 
sponds to the evolution of volume filling fraction of Helll 
regions due to a late period of Hell reionization, driven by 
the population of known quasars a,t z > 3. Middle panel. 
The electron scattering optical depth, integrated from to 
z. Two barely distinguishable curves are shown here: ex- 
cluding/including free electrons from the quasar-induced 
Hell ionization. Bottom, panel. The solid curve shows 
the probability distribution of hydrogen reionization red- 
shifts. The dashed curve shows the distribution of (late) 
Hell reionization redshifts. The amplitude of the latter 
distribution has been divided by a factor 20, for clarity of 
presentation. 



are impervious to photoionization feedback, and complete 
the reionization of hydrogen. 

As far as the late time thermal state is concerned, the 
relevant ionizing spectrum is that of the sources that turn 
on after 2 ~ 17. These are Population II stars - a rea- 
sonable spectrum is am = 1, and auei = 4, and heavily 
truncated beyond i^Heii (Leitherer et al. 1999). As we 
have illustrated with examples in the last section, a spec- 
trum as soft as this, even with some amount of late (HI, 
Hel but not Hell) reionization after z ~ 10, has difficulty 
heating up the IGM to high enough temperatures (even 
after accounting for processing of spectrum by the IGM, 
see §3). 

We are therefore led to consider the effect of quasars. 
Sokasian, Abel & Hernquist (2002) showed, based on a 
3D radiative transfer simulation, that Hell reionization 
should occur around 2; ~ 4, driven by the population of 
known quasars. ^ According to this work. Hell— >HeIII 

® Wyithe & Loeb (2002) considered Hell reionization by the known 
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Fig. 6. — Probability distribution of temperature To for 
fluid elements at mean density, for three different redshifts 
z — 2.4 (blue dashed histogram), 3.0 (magenta dotted his- 
togram), and 3.9 (red solid histogram). The observed tem- 
peratures with their 2 a errorbars are shown at the top. 



reionization completes within a relatively short redshift 
range (see the left dashed curve in the top panel of Fig. 
5). We approximate their probability distribution of (late) 
Hell reionization redshifts as a Gaussian centered at 4.1 
with a fuU-width-at-half-maximum of 0.5 (dashed curve in 
the bottom panel of Fig. 5). The quasar spectrum is as 
described in §3. 

We generate reionization histories for an ensemble of 
fluid elements (at mean density) in a stochastic way: de- 
termining the redshifts of HI (and Hel) reionization, and 
Hell reionization by drawing from the probability distri- 
butions described in the bottom panel of Fig. 5. The 
resulting temperatures at 2; = 2.4, 3.0, and 3.9 have a 
scatter because of the stochastic history, and their proba- 
bility distributions are shown in Fig. 6. 

The results shown in this figure are intriguing. At 
z = 3.9, the distribution (red solid histogram) peaks 
around To = 2.2 x 10** K, but it has a significant bump 
below 10** K as well. This bimodal distribution is a result 
of the stochastic reionization history: at z = 3.9, there 
is a minority of fluid elements that have not undergone 
Hell reionization, and so they are significantly colder, by 
more than a factor of 2. There have been attempts in the 
past to look for temperature fluctuations (beyond what 
one expects from variations with density) in the Lya forest 
(Zaldarriaga 2002, Theuns et al. 2002b). So far there has 
been no detection. Our results indicate that the stochas- 
ticity of reionization can detectably increase the width of 
the temperature distribution. It will be very interesting to 

quasar population in a semi-analytical model and obtained a similar 
result: Hell reionization extends beyond ^ ~ 4, but is complete by 
around then. 



confirm our results by a more detailed analysis that folds 
in a realistic density distribution, and to apply the obser- 
vational search techniques to larger datasets and to higher 
redshifts (z 3 so far). 

The formal mean and rms scatter of Tq at z = 3.9 is 
1.8 x lO'' ± 7 x 10^ K. It can be seen from Fig. 6 that 
the scatter gets progressively smaller as one goes to lower 
redshifts: at z 3.0, To = 1.7 x 10"* ± 8 x 10^ K; at 
z = 2.4, To = 1.5 X lO'^ ± 5 x 10^ K. The scatter at these 
lower redshifts are no larger than the expected scatter from 
shock-heating as well as dynamics (e.g. Croft et al. 1997, 
Hui & Gnedin 1997, Dave & Tripp 2001). 

6. DISCUSSION 

We find that the temperature of the IGM, as inferred 
from Lya forest spectra at redshifts z « 2 — 4, leads to 
several interesting conclusions. Especially interesting are 
the conclusions we obtain when the Lya forest tempera- 
ture is considered together with the WMAP results. Our 
main conclusions are summarized as follows: 

• A vanilla reionization model, where Xm, Xuci, and/or 
Xneii undergo order unity changes at z > 10, but 
suffer no such changes at z < 10, is ruled out by 
temperature measurements of ZHTOl, especially at 
z = 3.9. This relies on assumptions about the ion- 
izing spectrum, which are discussed in §3. For a 
power-law spectrum, rigorous requirements can be 
put on the hardness of the ionizing background to 
evade the the above argument: a < 0.12 if helium is 
doubly ionized, or a < —2.2 if helium is singly ion- 
ized, for J u~°' (with a cutoff at i^hcII if helium is 
only singly ionized) . The idea of an asymptotic evo- 
lution for the thermal state of the IGM is quite useful 
in formalizing the above argument. The asymptote 
(at 2 < z < 4) is described quite accurately (~ 5%) 
by eq. (4). It is reached as long as any order unity 
changes in Xhi,^hcI, and/or XhoII occur prior to 

z = 10. 

• Conversely, the requirement by WMAP that the uni- 
verse reionizes early, at z > 10, implies the reion- 
ization history is complex. To fulfill the temperature 
constraints, there must be additional periods of or- 
der unity changes in one or more of the fractions 
^Hii^Hcij ^Hoii O't redshift below 10. This is an 
argument separate from the argument based on the 
evolution of mean transmission (Haiman & Holder 
2003; based on the Gunn-Peterson troughs observed 
by Becker et al. 2001 and Fan et al. 2003 at z 6, 
and comparison against an extrapolation from lower 
redshifts). 

• Exactly what kind of changes at z < 10 are neces- 
sary to match the temperature constraint at z — 3.9 
depends critically on the ionizing spectrum. If the 
spectrum is harder than i^^^-^ above the HI and Hel 
ionization thresholds (i.e. harder than the 'quasar' 
spectrum discussed in §3), then order unity changes 
in Xm and Xuei in the period 6 < z < 10 (the lower 
limit of 6 being set by the SDSS mean transmission 
measurements; Fan et al. 2002), even without ac- 
companying changes in XhcII (i-e. no Hell reioniza- 
tion), are sufficient to heat up the IGM to within the 
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constraint at z — 3.9. Conversely, if the spectrum is 
not hard enough (for instance, a population II stel- 
lar spectrum described in §5), some amount of Hell 
reionization-heating is necessary prior to 2; = 3.9. 

• The IGM temperature can have a broad distribution 
(beyond what one expects based on variation luith 
density) close to reionization events, but the scat- 
ter diminishes with time. The reionization history 
of the universe is almost certainly stochastic, in the 
sense that different fluid elements reionize at dif- 
ferent times, depending on when a given element 
becomes engulfed in expanding HII (or Helll) re- 
gions around ionizing sources. Based on simple semi- 
analytic models, we predict the probability distribu- 
tion of ionization redshifts. The resulting tempera- 
ture scatter is particularly large close to reionization 
events, but diminishes quickly with time. This is 
illustrated in Fig. 6. It would be very interesting 
to search for the broad initial temperature scatter 
that is predicted by physically motivated reioniza- 
tion models. 

• The known quasar population at z < 5 may heat the 
IGM to sufficiently high temperatures by reionizing 
Hell. In §5, we describe a physically motivated reion- 
ization model with a stochastic history. Hydrogen 
(and Hel) reionization is accomplished early on by 
metal-free stars (17 <z< 30), and later on by pop- 
ulation II stars (z < 17), and completes by z ^ 6.5. 
This fulfills the dual requirements of a large CMB op- 
tical depth, and the Gunn-Pctcrson trough seen by 
SDSS. Because the population II stars have a rather 
soft spectrum, photo-heating of HI and Hel alone 
is not enough to bring the temperature up to the 
observed level, and this model fails to predict sufh- 
ciently high IGM temperatures - despite the percola- 
tion occurring at the relatively low redshift of 2; ~ 7. 
We therefore make use of the predictions of Sokasian 
ct al. (2002) for a brief period of Hell reionization 
around z ~ 4, based on the known population of 
quasars. We find that such a model can satisfy the 
temperature constraints. It is important to note, 
however, that this is not the only way to achieve the 
observed temperatures. For instance, a population 
of dim quasars can turn on at 2; > 6 (Haiman & 
Loeb 1997), which has a sufficiently hard spectrum 
to either heat the IGM via photo-heating of HI and 
Hel alone ( but their spectrum has to be harder than 

if Hell is not reionized), or boost the temper- 
ature via Hell reionization as well. If Hell reioniza- 
tion takes place via these mini-quasars, the turn-on 
of the known population of quasars at z < 4 then 
has a much weaker influence on the thermal state at 
low redshifts. 

• Hell reionization makes little difference to the elec- 
tron scattering optical depth, particularly if it occurs 
late. For instance, in the model shown in Fig. 5, 
extra electrons from Hell reionization changes Tg at 
a 2% level fractionally. 



possibilities are galactic outflows and Compton heating 
by a hard X-ray background. Adclbcrgcr ct al. (2002) re- 
cently observed signatures of galactic winds into the IGM. 
Such outflows can in principle heat up the IGM. However, 
observations by Ranch et al. (2001) of close pairs of lines 
of sight in lens systems suggest that the IGM is not turbu- 
lent on small-scales, arguing against significant stir-up of 
the IGM by winds. Moreover, galactic outflows can heat 
up the IGM to a variety of temperatures - the fact that 
the observed temperatures are in the range of expecta- 
tions for a photo-ionized, and photo-heated, gas suggests 
photo-heating is the simplest explanation. 

Another important question is whether Compton heat- 
ing by a hard X-ray background (XRB) could be more 
important than photoelectric heating. This question was 
considered by Madau & Efstathiou (1999), who assumed 
that the redshift evolution of the sources of the hard XRB 
parallels the flat distribution that had been determined for 
the soft X-ray AGN luminosity function beyond 2; ~ 2. 
Under this assumption, they found the hard XRB to be 
an important source of heating relative to the UV back- 
ground at redshifts 2 > 2, raising the IGM temperature to 
1.5 X lO'^K at 2 ^ 4. Two new developments make it un- 
likely for the hard XRB to be an important source of heat- 
ing. First, as pointed out by Abel & Haehnelt (1999), the 
photoelectric heating rate can be signiflcantly increased in 
optically thick gas (when hydrogen and helium flrst gets 
ionized); we here adopt these increased rates. Second, the 
sources of the hard XRB have been resolved by the Chan- 
dra satellite, and, rather than paralleling the soft X-ray 
luminosity function, the hard X-ray sources exhibit a steep 
decline towards high redshift beyond 2 ~ 2 (Cowie et al. 
2003). 

Our investigation raises a number of interesting issues. 
Is Hell reionization prior to redshift 4 necessary to heat 
up the IGM to the right level? This depends critically on 
what kind of sources and spectra are available. The pre- 
diction for a large scatter in temperature (factor of about 
2) close to the epoch of Hell reionization is something 
one could look for, if indeed Hell reionization occurs late 
(Zaldarriaga 2002, Theuns et al. 2002b). In this paper, we 
have focused entirely on the thermal state of fluid elements 
at the mean density. The formalism of Hui & Gnedin 
(1997) can be used to compute the same for elements at 
Sp/p < 10. It is interesting to explore how the variation of 
temperature with density (T oc p'^^^ . an effective equation 
of state) might place additional constraints on the reion- 
ization history. At present, measurements of 7 are quite 
noisy (e.g. ZHTOl). New approaches to constrain it better 
will therefore be very useful (Dijkstra, Lidz & Hui 2003). 
Finally, there is the issue of exotic ionization mcchanisins 
(e.g. Berezhiani & Khlopov 1990) which is particularly 
interesting in the face of early reionization. 

We thank Renyue Ccn, Andrea Fcrrara, Wayne Hu, Avi 
Loeb, Jerry Ostrikcr, Joop Schaye, and especially Adam 
Lidz, for useful discussions. LH is supported in part by an 
Outstanding Junior Investigator Award from the DOE, 
an AST-0098437 grant from the NSF, and by the DOE at 
Fermilab, and NASA grant NAG5-10842. 



An important caveat in our discussion so far is the possi- 
ble existence of additional heating mechanisms. Two such 
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